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Abstract 



. . . On the basis of systematic measurements of fragmentation reactions, which provide a detailed overview on the velocity 
distributions of residual nuclei, an improved description of the kinematical properties of the fragmentation residues 
is established. This work is dedicated to the fluctuations of their momentum distributions. In contrast to previous 
investigations, limited to close-to-projectile fragments, we extended our study to the entire production range, down 
"to the lightest observed fragments. In this context, beside the contribution of abrasion and evaporation processes, 
■we considered the effect of the thermal break-up on the width of the momentum distributions. Using approximated 
"theoretical descriptions of the different reaction stages, a new analytical formula for the variance of the momentum 
distribution is derived, which is well adapted to technical applications. 
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1. Introduction 

Heavy-ions reactions with various target materials at 
energies well above the Fermi-energy regime, have gained 
■increasing interest for a variety of applications. They are 
used for the production of rare isotopes in existing [iHll 
and future secondary-beam facilities [B43> they are ex- 
ploited in the field of nuclear technology for the opera- 
tion of spallation neutron sources in particular for 
the incineration of nuclear waste in an accelerator-driven 
system they are a source of activation due to in- 
evitable beam losses along high-energy accelerators and 
•in the beam dump ^llj . Besides their cross sections, the 
kinematical characteristics of the fragmentation residues 
emerging from these reactions are very important. In in- 
flight secondary-beam facilities, they determine the emit- 
tance of the rare-isotope beams; in irradiated construction 
material, they are responsible for aging phenomena. 
1 In non-central collisions of heavy ions at energies far be- 
yond the Fermi-energy regime, the geometrical abrasion 
model fl^l suggests a clear cut between participant mat- 
ter, corresponding to the overlap zone of the two nuclei, 
and spectator matter outside the overlap zone. While the 
participant matter is subject to strong heating and com- 
pression due to high-energy nucleon-nucleon collisions, the 
spectator matter is sheared off from the projectile, respec- 
tively target nucleus, and continues moving on with es- 
sentially its initial velocity. On a smaller scale, however. 
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a shift in the mean longitudinal velocity of the spectator 
matter is observed. For very peripheral collisions, fric- 
tion dominates , while for less peripheral collisions an 
opposite effect has been observed [14|, which has been at- 
tributed to the response of the spectator to the participant 
blast In addition, the velocity distribution of the 

spectator matter is characterized by fluctuations around 
its mean value. The magnitude of these fluctuations is 
most often larger than the shift of the mean value and, 
thus, constitutes the most important characteristics of the 
kinematical properties of the spectators. 

The basic theoretical understanding of the fluctuations 
of the momentum distributions dates back to the pioneer- 
ing work of Goldhaber in 1974 [l^l, who described two 
possible scenarios, the resulting recoil momentum due to 
a sequential particle emission, or the random superposition 
of the individual Fermi momenta of the nucleons removed 
in an abrasion process. A quantitative simulation of the 
velocity distributions of individual nuclides, also concen- 
trating on very peripheral collisions, considering the con- 
tributions of the abrasion and the evaporation stages by a 
full Monte-Carlo nuclear-reaction code, has been presented 
by Hanelt et al. • They have shown that the recoil mo- 
mentum of evaporated particles modifies the momentum 
dispersion of final fragments. A comprehensive system- 
atics on momentum distributions from reactions with rel- 
ativistic ions, documenting the status of that time, was 
established by Morrissey in 1989 |18i], and an empirical 
formula for very peripheral collisions was proposed: 
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where Ap and A f represent the mass of the projectile and 
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the mass of the final fragment, respectively. Although this 
empirical formula perfectly reproduces most data near the 
projectile, it becomes unrealistic for lighter fragments. In- 
deed, the predicted dispersion prediction does not decrease 
towards low masses and finally gives unrealistically large 
values for very light fragments. One has to realize that 
the width of the momentum has to go down to zero with 
mass approaching zero! Thus, this formula cannot cover 
the whole range of fragments produced and observed in 
fragmentation reactions. 

Recently, a semi-empirical model describing the momen- 
tum distribution of the fragmentation residues has been in- 
troduced by Tarasov ^19J . The momentum distributions is 
represented by a convolution of a Gaussian with an expo- 
nential tail, where the shape of the convoluted distribution 
and its variation with projectile energy and the masses of 
projectile and fragment are associated with the influence 
of abrasion, friction and evaporation on the kinematics of 
the fragments. The free parameters of that distribution 
are tuned to fit experimental data. 

The need to build up a model fully based on theoret- 
ical considerations, taking into account the multiple pro- 
cesses at play in fragmentation reactions, valid over the 
whole mass range of produced residues but compact and 
fast enough to be used for various applications, is then 
obvious. 

Thanks to a series of extensive experiments on system- 
atic measurements of fragmentation residues, mostly per- 
formed at GSI, Darmstadt, the empirical knowledge in this 
field has improved substantially [17 , By investigat- 



ing the projectile-like fragments, a detailed overview on 
cross sections and velocity distributions of individual nu- 
clides has been obtained. The high quality of these data 
comes not only from the high precision, but also from the 
fact that, in contrast to experiments on the properties of 
target-like fragments, which suffer from a low-energy cut- 
off, the kinematical properties of the projectile-like frag- 
ments can be fully measured. In addition, for a few reac- 
tions, the final residues were measured over a very broad 
mass range. The availability of new, high-quality data has 
motivated us to re-visit the understanding of the kinemat- 
ical properties of the fragmentation residues. The present 
work is dedicated to the fluctuations of their momentum 
distributions. 

It is our aim to study the fluctuations of the momentum 
distributions of all fragmentation residues emerging from 
the spectator matter. In contrast to many previous inves- 
tigations, we consider the whole mass range, from the pro- 
jectile, respectively target nucleus, down to light charged 
particles. 

For this purpose, we also treat the influence of the multi- 
fragmentation stage on the kinematical properties of frag- 
mentation residues. This break-up of the system due to 
thermal instabilities becomes important at high beam en- 
ergies and was not considered in previous works. 

In the present paper, we develop a model fully based 
on theory. We derive analytical approximations from ap- 



propriate theoretical descriptions of the different stages of 
the collision. This is used to develop a new comprehen- 
sive analytical formula for the variance of the momentum 
distributions, well adapted to technical applications. 



2. Experimental data 

From the body of data on kinematical properties of frag- 
mentation residues, we have chosen three experiments per- 
formed at the FRagment Separator (FRS) in GSI, Darm- 
stadt. They represent the fragmentation of two light nu- 
clei (^^Fe [H3|, ®^Kr and a medium-heavy nucleus 
(^■^^Xe [H^). These experiments cover a broad mass range: 
residues have been observed from the vicinity of the pro- 
jectile down to elements lighter than neon. 

The FRS is characterized by a high momentum reso- 
lution and a limited acceptance both in momentum and 
angle [gH. The momentum resolution (FWHM) amounts 
to Ap/p ~ 5 • 10^^, the momentum acceptance is ±1.5%, 
while the angular acceptance comprises 15 mrad around 
the beam axis. The limited momentum acceptance is not 
crucial; since the full momentum range can be covered 
by combining results obtained for different settings of the 
magnetic fields. We limit the present work to the lon- 
gitudinal momentum distribution of the residues emitted 
within the angular acceptance of the FRS. Nevertheless, 
the same formalism can be used for the transverse compo- 
nent of momentum distributions |2C 



62| . Since at the pro- 



jectile energies considered in the present work, the produc- 
tion of fragmentation residues is approximately isotropic 
(20l . [b^ . [63| , our description should be valid for the trans- 
verse momentum dispersion as well. 

Three measured longitudinal velocity distributions, in 
the projectile frame, are given as examples (fig. [1]). They 
correspond to ^^''Sn, ^°Ni and ^'^C fragments from the re- 
action ^■^^Xe+Pb at 1 ^ GeV [56,]. Over a large mass range, 
the residues follow a Gaussian distribution. Only the light- 
est residues show a more complex distribution with an 
asymmetry or even several humps. Light nuclei can origi- 
nate from two different processes (iol . 53, 64, 6^. In binary 
mass splits, e.g. very asymmetric fission. Coulomb forces 
lead to a two-humped longitudinal velocity distribution 
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49|. In case of multifragmentation, the Gaussian shape 
of longitudinal velocity spectra is preserved. Please note 
that our model aims at describing fragmentation reactions, 
and therefore cannot be used to calculate the momentum 
width of fission residues. 

The data from ref. (23 . [s^ , [56| will be compared with 
the predictions of our analytical model in the section HI 
but let us first present its fundamental ideas. 



3. Pragmentation mechanism 

Let us remind the different descriptions of the processes 
that happen during the collision and the de-excitation 



2 



3.1. Abrasion 




Figure 1: Longitudinal velocity distributions in the projectile frame 
for various fragments (transmitted through the 15 mrad angular ac- 
ceptance of the FRS) from the reaction ^^®Xe+Pb at 1 GeV per 
nucleon, taken from (56| . Heavy fragments follow a Gaussian distri- 
bution. On the contrary, the lightest ones populate an asymmetric 
distribution. 



phase of the fragmentation reaction and study their in- 
fluence on the width of the longitudinal momentum distri- 
butions of the fragments. 

As a typical collision configuration, we will assume a 
projectile nucleus hitting a target nucleus (at rest) with 
an impact parameter in principle different from zero. In 
the overlap region, projectile and target nucleons interact 
strongly, while the rest of the projectile and the target nu- 
cleus, so called spectators, are almost undisturbed by the 
collision. This picture is well suited in a wide energy do- 
main above the Fermi energy, and the model developed in 
the present work can thus be used for comparison with var- 
ious projectile- and target-fragmentation data from avail- 
able devices (e.g. at GSI, MSU, RIKEN. . . ). In the fohow- 
ing, we will concentrate on the projectile-like fragments, 
since we will compare the predictions of the present model 
with the projectile-like residues measured at the FRS in 
GSI. Of course, the model is also applicable to target-like 
residues. We will discuss the contributions of the differ- 
ent stages of the fragmentation process to the momentum 
dispersion of the final residues. 



Dedicated experiments have long shown that the lon- 
gitudinal momentum of the heavy fragmentation residues 
follows a Gaussian distribution in the projectile frame. A 
few models have been developed to describe the standard 
deviation of this distribution. 

For the standard deviation of momentum in longitudi- 
nal or transverse directions {(^p^^ and <Jpj,), the most of- 
ten used theoretical model is Goldhaber's prediction for 
the dispersion induced in the abrasion process due to the 
Fermi motion of the nucleons [16]. 

A proposition of Goldhaber is to consider the reaction 
as a sudden cut-off of a part of the projectile, without 
taking into account any further evolution of the remain- 
ing part of the projectile. Some nucleons are removed in- 
stantaneously, without inducing momentum transfer. This 
reaction mechanism is referred to as abrasion, and it sug- 
gests already a removal of matter by friction phenomena, 
so that its description in [f 6] is probably too simple, but we 
shall discuss that after considering the strong implications 
of the assumption made by Goldhaber. An instantaneous 
removal of several nucleons from the projectile does in- 
deed affect the dynamical features of this projectile. In 
fact, this contribution explains great part of the measured 
momentum width of the fragmentation residues. The re- 
moved nucleons are well defined by their positions inside 
the nucleus at the moment of collision. Yet, the momenta 
of these nucleons are sampled over a broad distribution. 

Considering the ensemble of nucleons as a Fermi gas, 
we know that these constituents of the nucleus have an in- 
trinsic movement, even at zero temperature. This internal 
motion affects the observable features of the projectile-like 
fragments. In a Fermi gas inside a square- well potential 
with infinitely high borders, there is no correlation be- 
tween the position and the momentum of a nucleon. Thus, 
the momentum of an abraded nucleon is a random sample 
from the projectile's Fermi sphere filled up to the radius 
Pf, in Goldhaber's view [31 . The authors of the present 
paper are aware of the works of Friedman (66| demon- 
strating that absorption, in peripheral collisions, prevents 
from sampling over the whole nucleus. Nevertheless, as 
one can see from ref. [ill, the predictions with the absorp- 
tion taken into account are rather close to the Goldhaber 
model. Therefore, and owing to the transparency of that 
latter, it was decided to use Goldhaber model as a basis 
of the present work. 

The mean square momentum of a nucleon in a Fermi 
gas is 
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The momentum variance ap of the projectile spectator 
is given by the sum of the individual contributions from 
the abraded nucleons. Provided the assumption that this 
broadening is equally distributed in all directions, the pro- 
jection along the beam axis gives — (Tp/3 (it is also 
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relevant for transverse momentum studies, see for exam- 
ple IHi). 

The variance is then linked with the Fermi momentum 
Pf of the projectile of mass Ap and the mass of the frag- 
ment after the abrasion process, so-called prefragment, A: 



p| AjAp - A) 
5 ■ (Ap-l) 



(3) 



The Fermi momentum varies only weakly with the mass 
of the nucleus. One can find an expression of the mass 
dependence of the Fermi momentum in reference |67| that 
relies on data from 16811: 



PFiA) = 281 • (1 - yl-"-^*'**) MeV/c . 



(4) 



We can comment that eq. [3] is symmetric regarding the 
mass, so that an abrasion process that removes half of the 
nucleons in the projectile produces the largest width in 
the distribution. In this hypothesis, light fragments and 
fragments of mass close to that of the projectile should 
have a similar broadening of their longitudinal momentum 
distributions. 

This model reproduces some data rather well. At the 
time of release of Goldhaber's article, most data available 
were limited to small mass losses, corresponding to the 
heaviest fragments. Goldhaber's formula is indeed in quite 
good agreement with these data. 

But Goldhaber's model is incomplete in the sense that 
the spectator matter after the abrasion process is excited. 
Therefore, this model does not refer to the measured frag- 
mentation residues, which result from a subsequent deex- 
citation stage, but rather to an excited prefragment. The 
comprehension of the deexcitation processes is the key to 
address a model for the final fragment. 

3. 2. Evaporation 

The removal of nucleons by abrasion is a step that does 
not leave the nucleus cold. It can be seen as a friction 
process, removing nucleons from the projectile, but also 
inducing excitation energy in the remaining fragment. It 
was shown by different methods (23l . [g^ that this energy 
induction amounts on average to 27 MeV per abraded nu- 
cleon. Therefore, the spectators may be highly excited, 
giving rise to the evaporation of a considerable number of 
nucleons and light nuclei. 

The evaporation stage can be described by an appropri- 
ate code in the frame of the statistical abrasion-ablation 
model pTOl, using realistic binding energies. Coulomb bar- 
riers and level densities. Each emitted particle induces a 
recoil momentum to the corresponding residue according 
to its kinetic energy and momentum conservation. The in- 
fluence of the evaporation process has been demonstrated 
by Hanelt et al. [17], stating that the recoil-momentum 
induced by sequentially evaporated nucleons should no- 
ticeably modify the width of momentum distributions. 

In most cases, the individual contributions of the evap- 
orated particles are small and just slightly increase the 



width of the Gaussian distribution resulting from the abra- 
sion stage. Only in specific cases, like symmetric fission, 
the recoil momentum is so large that it dominates the 
kinematic properties of the final fragments, and devia- 
tions from a Gaussian distribution are observed. Although 
fission plays an important role for very heavy fragments, 
we will consider only non-fissioning systems in the present 
work. 

De Jong et al. [7l[ included the influence of the evapo- 
ration process on the angular momentum distribution in 
an analytical formula. The basic idea was to assume a 
certain proportionality between abraded mass and evapo- 
rated mass. They obtained satisfactory predictions with a 
single parameter. They used 77 = 2 as a typical value of 
mass loss due to evaporation per abraded nucleon. This 
means that the abrasion of one nucleon leads on average 
to a mass loss of two nucleons by evaporation. Their de- 
scription was formulated as: 
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We have adapted this idea for expressing the conse- 
quence of evaporation on the linear momentum by a dedi- 
cated analytical expression. It will be presented in section 
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For the prefragment, the release of nucleons (mostly in- 
dividually but also as bound clusters) decreases its exci- 
tation energy. Taking into account the mean separation 
energy of one nucleon and the average nuclear tempera- 
ture inside the prefragment, this "cooling" of the system 
allowed by evaporation is considered to be 15 MeV per 
emitted nucleon on average. 

Since the angular momentum induced in an abrasion re- 
action is rather low [tH , evaporation from the excited pre- 
fragment is an isotropic process in the frame of the source 
of emission. Thus, the mean velocity of the prefragment is 
not affected by evaporation. The momentum, however, is 
reduced on average because of the mass loss. In addition, 
the recoil momenta of the individual evaporated particles 
tend to increase the width of the momentum distribution. 

In order to compare explicitly the impact of the recoil 
process with the Fermi contribution discussed before, let 
us deflne , the ratio between the mean individual en- 
ergy taken from the prefragment by the evaporation of 
one nucleon < -Ecvap > and the mean individual energy of 
one nucleon inside the nucleus due to the Fermi motion 
< Eu 
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The mean kinetic energy < -Efcrmion > of a nucleon in a 
Fermi gas is: 
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The average Fermi energy Ep of a nucleus amounts to 
about 33 MeV; this value leads to < i?formion >— 20MeV. 
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The value of < -Ecvap > is not easily estimated, because 
it depends on the charge of the mother nucleus, and of 
course is different for each type of evaporated particle (neu- 
tron, proton, alpha, intermediate-mass fragments...). For 
simplicity, we consider only the evaporation of neutrons 
and protons. Moreover, we assume that they are emit- 
ted with equal probability, so that the mean recoil energy 
< £^cvap > is given by the mean value between neutron and 
proton kinetic energies En and Ep , respectively. We fix the 
thermal contribution to En and Ep to 8 MeV; for protons, 
we estimate the contribution from Coulomb repulsion us- 
ing the mass Ap and nuclear charge Zp of the projectile. 
In this way, we obtain for 77^ the following expression: 
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with ro=1.4 fm and e the elementary charge. 

The mean square momenta of a single nucleon, for either 
Fermi-gas contribution or for the evaporation recoil, are 
directly proportional to the corresponding mean energy: 



< p >= 2m„- < E > 



(9) 



with mn the mass of a nucleon. The ratio between the 
mean energies is then equivalent to the ratio of momenta 
squared: 
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using the expression of < p|,j.,iiion ^ from eq. ^ For the 
mean square momentum in the longitudinal direction, one 
finally has: 



< P||cvap > = 



(11) 



We shall use this expression to take the recoil momentum 



into account in the formula for ct^ (see details in section 



"analytical model") . 

3.3. Multifragmentation 

It is commonly accepted that a nucleus starts to show 
some thermal instabilities when its excitation energy ex- 
ceeds a value of about 3 MeV per nucleon. This value 
can easily be reached through the abrasion process. The 
kinematics of the fragments is defined when they cease 
to interact, at freeze-out. The momenta of the fragments 
are given by a radial expansion, a random motion of the 
nascent fragments and by a consecutive Coulomb expan- 
sion. 

3.3.1. Random motion 

For estimating the random motion, two rather diverg- 
ing ideas were introduced in the literature. In one extreme, 
it was assumed that the fragments behave like molecules 
of an ideal gas [72|. In this case, their individual mo- 
tion is governed by a Maxwell-Boltzmann distribution and 



the variance of the momentum distribution of a fragment 
of mass A produced in the break-up of a system of mass 
Abu by multifragmentation is given by the following for- 
mula [161] : 



CT,, = mnkeTbu 



A(Abu - A) 



A 



bu 



(12) 



where Tbu is the freeze-out temperature, ke is the Boltz- 
mann constant and mn denotes the nucleon mass. 

In the other extreme, the nascent fragments keep the 
Fermi motion of the individual nucleons in the common 
source [sHj. In this second case, the distribution follows 
the same functional form, but the temperature parameter 
cannot be interpreted as a real temperature. 

Odeh et al. have shown that the apparent tempera- 
ture Tkpp extracted from the kinetic-energy spectra of the 
multifragmentation products is much larger than the tem- 
perature extracted using isotopic or excited-states popula- 
tion thermometer [35j . Such large apparent temperatures 
would indicate that the process of simultaneous break-up 
is rather fast and that the system does not have time to 
reach full thermal equilibrium. In other words, the Fermi 
motion inside the multifragmenting source is mostly re- 
sponsible for the broadening of the momentum distribu- 
tion of the multifragmentation products. Thus, ideas be- 
hind the Goldhaber formula (eq. [3]) should also be valid 
for describing the momentum dispersion of the multifrag- 
mentation products. However, one should not forget that 
in case of multifragmentation the multifragmenting source 
is not cold as in case of abrasion, but has a finite temper- 
ature Tbu- Thus, the mean velocity of the nucleons inside 
the multifragmenting source is larger than the mean veloc- 
ity of the nucleons inside the cold source. Consequently, 
the slope parameter (i.e. apparent temperature) of the 
kinetic-energy spectra of the multifragmentation products 
is increased relative to the case of the cold fragmenting 
source. This effect has been studied by W. Bauer [73], 
and he developed an analytical expression for calculating 
the momentum dispersion of a multifragmenting Fermi gas 
at finite temperature. According to ref. [73| . the apparent 
temperature Tapp can be expressed as a function of the 
real source temperature Tbu through a Taylor expansion 
by: 
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with Abu the mass of the multifragmenting source, Ep its 
Fermi energy, and A the mass of a produced fragment, 
which will enter the evaporation stage. This formula is 
incorporated in our model. 

3.3.2. Radial expansion 

The effect of the radial expansion is a lowering of the 
Fermi energy Ep and of the Fermi momentum pp down 
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to values -Epbu and PFbu according to the following expres- 
sions: 



Ephu — Ep ■ ( — ^ 



PFbu = Pf ■ 
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(15) 



In above equations, Vq is the normal volume, while Vbu 
is the volume at freeze-out. The generally accepted vol- 
ume increase at freeze-out is still not well established, and 
values between one and six times the normal volume are 
often considered [t^I- 

3.3.3. Coulomb expansion 

Furthermore, we shall consider the fact that the clus- 
ters formed in the thermal break-up are charged particles. 
Thus, Coulomb forces between these fragments have to be 
taken into account. 

K. C. Chung et al. developed ideas regarding the effects 
of this Coulomb expansion in nuclear fragmentation [74]. 
We summarize here the main ideas and results of their 
article. 

Assuming that the expansion is uniform, but that the 
fragment radii do not change during expansion, it is possi- 
ble to express the Coulomb contribution to the final kinetic 
energy: for a fragment of nuclear charge Z and mass A, 
at a distance r from the center of a prefragment of mass 
^prc, nuclear charge Zpio and radius i?pro, one has: 
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Equation [in] is remarkably simple, especially owing to the 
fact that the only relevant initial condition is the relative 
position of the fragment {A,Z) inside the system at break- 
up. 

While that kind of contribution can easily be included 
in a Monte-Carlo simulation, its inclusion in an analytical 
model may be more complex, or require some assumptions 
that we will discuss in the next sections. 

3.4-. Overview on the different fragmentation stages 

The impact parameter of the collision determines the 
mass and the excitation energy of the prefragment and 
hence the mass of the final fragment [t^. Therefore, one 
can make the connection between the mass of the final 
residue and the conditions at the beginning of the evapo- 
ration stage [s^l- This is schematically shown in fig. [21 

The different stages of the reaction are connected with 
a reduction of mass. To follow the evolution of a prefrag- 
ment, one has to read fig. [5] from right to left. The first 
process of the fragmentation reaction, the abrasion stage, 
reduces the mass of the projectile, respectively target nu- 
cleus, inducing substantial excitation energy. The average 
excitation energy follows the line which corresponds 



breakup 



abrasion 

evaporation 




Figure 2: Evolution of the excitation energy of the fragment through 
different reaction stages, from a given projectile of mass Ap. The 
abrasion process is represented by the line iJJ* = A_E ■ {Ap — A). If 
multifragmentation occurs, the system after freeze-out lands on the 
line E2 = A/K ■ T^^. In all cases, evaporation goes on average along 
■ [A — Af). See text for details. 



to AE=27 MeV induced by abrasion of one nucleon, as 
mentioned above. 

If the spectator after abrasion is left with an excitation 
energy of more than about 3 A MeV, it breaks up into 
several fragments simultaneously. The products of this 
multifragmentation process are formed with a freeze-out 
temperature Tbu of about 5 MeV (s^l- Before entering the 
evaporation process, these fragments fall on the freeze-out 
line £'2 ■ The relation between temperature and excitation 
energy is described by the following formula: 



E* = aT^ 



(17) 



In this equation, a is the level-density parameter; it is 
evaluated in a simplified way as A/K MeV~^. With a 
value of the inverse level-density parameter K=ll MeV 
and a freeze-out temperature of 5 MeV, we obtain E2 = 
ff ^ MeV. 

As soon as the system is excited (after abrasion or after 
multifragmentation) above the particle emission threshold, 
it has the possibility to emit neutrons, light charged par- 
ticles and light nuclei by evaporation. This deexcitation 
process is represented by the arrows . The link between 
the observed mass A f and the mass of the prefragment be- 
fore evaporation is given by E^ = e ■ (A — Af)^ where e is 
the parameter representing the average energy consumed 
per evaporated nucleon. We assume a value of e=15 MeV 
(see section [XU . 

Figure [5| suggests that final fragments lighter than a 
certain limiting mass An^ are products of multifragmen- 
tation, while fragments heavier than this limit should be 
regarded as simple abrasion-evaporation residues. 

In other words, the present model considers two regions 
in mass, corresponding to the two regimes preceeding evap- 
oration. They are divided by Aii,„ which is obtained by the 
projection following an evaporation line E^ of the intersec- 
tion of lines E* and E2 . The solution of this geometrical 
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requirement gives: 
K -e- 



rp2 



K -AE 



K ■ 



{K ■ AE- 



(18) 



Again, final residues with masses above An^n are the 
products of a pure abrasion-evaporation process, while fi- 
nal residues with masses below this limit have experienced 
a multifragmentation process. 

4. Analytical model 

Our goal is to obtain an analytical description of the 
momentum dispersion of fragmentation residues. Such an 
approach requires some approximations. 

We have seen that formula[3]represents only the abrasion 
stage of the fragmentation reaction and does not include 
the evaporation nor the multifragmentation stages. 

Nevertheless, formula |3] serves as a basis for the present 
work, since the Fermi motion is the dominant effect. 

Most of our results come from geometrical deductions 
of the reaction scheme depicted in fig. [21 We already ex- 
pressed the mass limit under which an observed fragment 
has presumably undergone break-up (eg. I18p . 

We denote the mass of the prefragment just before 
evaporation by Aq^ (GH stands for Goldhaber), since 
we will use it as a variable in Goldhaber's formula. In 
the abrasion-evaporation scenario, the mass before evap- 
oration is obtained at the intersection between lines Ei 
and £'1 (see fig. [5]). In the abrasion- multifragmentation- 
evaporation scenario, the mass before evaporation is the 
solution of the intersection of lines £'2 and E^ . For final- 
fragment masses close to the projectile {Af > A\i„^), we 
obtain: 



A 



GH 



eAf + AEAp 
{AE + e) 



and if Af < Ay, 



A 



K-e 



GH 



K-e 



T2 



■A, 



(19) 



(20) 



In the abrasion stage, the influence of the Fermi mo- 
tion is described by the Goldhaber formula (eq. [3]), and 
the standard deviation of the longitudinal momentum dis- 
tribution of abrasion-evaporation products including the 
change in mass due to evaporation is given by the follow- 
ing formula: 



PlIFo, 



A 



GH 



^gh(^p -^gh) 
{Ap - 1) 



(21) 



The Fermi momentum pp in the above equation is di- 
rectly obtained from eqU) On the other hand, in the case 
of multifragmentation, one must consider that the Gold- 
haber formula (eq. ^ does not account for the volume 
expansion during the break-up nor for the finite tempera- 
ture of the multifragmenting source. The volume expan- 
sion during the break-up is described via equations 1 141 and 



I15L where we assume a freeze-out volume equal 1.5 times 
the normal volume. The influence of flnite temperature 
of the multifragmenting source is described according to 
Bauer's analytical model presented by eq. 1131 where we 
assume a value of 5 MeV for the temperature at freeze-out 

Since the descriptions of abrasion and break-up are both 
based on the Fermi momentum, we can adapt the expres- 
sion of PFbu, into an apparent Fermi momentum: 



PFbu = Pf ■ 



14u 



1/3 



57r2 / Tbu 



12 \i?Fbu 



(22) 



The fact that we use an apparent Fermi momentum in 
case of multifragmentation requires that we write the cor- 
responding contributions to the momentum dispersion in 
two terms: one contribution from the abrasion process 
(formation of a prefragment of mass Abu out of a pro- 
jectile Ap), and another contribution from the multifrag- 
mentation (formation of a product of mass Aqu out of an 
ensemble of Abu nucleons). 

In case of multifragmentation, Aqh is expressed by 
eq. 1201 but we do not know a priori the mass Abu of the 
source of multifragmentation residues. As a reasonnable 
estimate, one can take the heaviest fragment that can en- 
ter multifragmentation. Let us recall that the abrasion 
stage is inducing an increase of excitation energy in the 
prefragment, by the amount of about AE = 27 MeV per 
abraded nucleon. Once the critical point leading to multi- 
fragmentation is reached, further abrasion of a few nucle- 
ons should lead to a break-up into smaller pieces, tending 
rapidly to vaporization. Final multifragmentation residues 
larger than alpha particles are then most probably pro- 
duced by a prefragment of mass close to the theoretical 
maximum. 

We calculated the mass Abu of the heaviest system un- 
dergoing multifragmentation. It corresponds to the colli- 
sions in which the abrasion has introduced enough energy 
to excite the nucleus to the critical temperature of 5 MeV, 
and is obtained from the intersection of the lines E^ and 
E2 shown in flgure[2j 



Abu = 



K ■ AE 



K-AE + T^^ 



■Ar. 



(23) 



Thus, the total contribution from the Fermi motion in 
the abrasion and the break-up on the kinematics of multi- 
fragmentation products is written: 



a: 



P II Fermi 



Abu 
Agh 



2 2 
Pf 



Abu(Ap — Abu) 



2 ~2 
Pphu 



[Ap 1) 

^GH(Abu 



A, 



GHJ 



(Abu - 1) 



(24) 



with pFhu the apparent Fermi momentum at break-up, 
given by eq. (55] 

Apart from the effects of the Fermi motion, we men- 
tioned that in case of multifragmentation, the expansion 
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due to Coulomb forces should be implemented in the 
model. Let us reconsider eq. [Tni Many variables are 
subject to fluctuations from one event to another, but 
through simple assumptions, one can sim plif y this equa- 
tion. Isospin thermometer experiments |23l.l39l| have shown 
that the nuclear composition of prefragments just before 
the sequential evaporation steps is compatible with a sce- 
nario in which the prefragment keeps "in memory" the ratio 
Zp/Ap of the projectile. Both processes, abrasion and mul- 
tifragmentation, preserve this value on average. It means 
that both the mother nucleus and its considered daugh- 
ter have on average the same proportion of neutrons and 
protons. We express then the atomic number of these two 
nuclei in the following form: 



An 



(25) 



where A, Z design the features of each considered part of 
the nuclear ensemble. The effect of neutron distillation [tgI . 
77| on the mean N/Z ratio of the fragments, which is small 
compared to the effect of evaporation, was neglected. Let 
us stress that this ratio is not preserved by the evaporation 
process, so that this formula does not apply to the final 
fragment. 

The system entering multifragmentation was taken to 
be the heaviest, with its mass Abu given by eq. [231 its 
charge Zbu is determined by the fixed ratio Z/A. This 
heaviest mother nucleus will generate the biggest electric 
field. This assumption should give us an upper estimate 
of the Coulomb repulsion, and at least allow us to check 
whether the magnitude of such an effect is compatible with 
the observed velocity distributions. 

The distance r between the fragments at their forma- 
tion is also crucial. We assume that the production of a 
fragment has a uniform probability inside the volume of 
the system. The variance of such a distribution of proba- 
bility in one direction (in our case, along the longitudinal 
axis) is then the one of a uniform distribution between 
an interval [-rmax;+''max]; such variance takes the form 
(2ri„ax)^/12. The integration of the electrical field in the 
spherical charged source up to the bounds of this interval 
leads to a Coulomb contribution to the longitudinal mo- 
mentum that we name Pcoui- The Coulomb contribution 
to the variance of the spectra is then written: 



Coul 



(26) 



with the contribution to the momentum (see eq. I16p taken 
as: 



P 2 

— — • r . 



•(l-4^f(27) 

^bu 



with u = 931.5 MeV/c2 and = 1.44 MeV-fm. r^ax is 
the maximum distance between the centre of mass of the 
fragment and the one of the whole system, so it is the 
radius of the mother nucleus, reduced by the radius of the 



produced piece: 



— ^bu — ^GH- 



The further refinement in the description of the pro- 
cesses undergone by the observed fragment is related to the 
dynamics acting in the evaporation step. The sequential 
aspect of the evaporation gives rise to recoil momentum, 
broadening the momentum distributions. 

For this approach, we look first at the variance of the 
longitudinal velocity distributions. If we start from the 
prefragment, we may assume that the variance after evap- 
oration of one nucleon is modified by the induced recoil 
momentum; the additional contribution to the variance of 
the velocity distribution is: 



P. 



cvap 



(Agh - 1)2 ■ 



(28) 



We suppose that a single particle is evaporated with the 
mean momentum Povap- Then for a number n of emitted 
particles, we can write: 



2 . \ " 

'^cvap ' / J 

4=0 



1 



(A 



GH 



(29) 



For a large number Aqh (compared to the number of evap- 
orated particles n), the difference between two steps of the 
sum is rather small. The sum can then be approximated 
by an integral. With the upper bound of the integral being 
n = Aqu — Af, we obtain this integrated relation for the 
final fragment: 



Pc 



1 



1 



A 



GH 



Since for the momentum variance, we have: 



A 2^2 



(30) 



(31) 



we can deduce the following, expressing p^vap along the 
longitudinal axis as in eq. [TT] (with 77^ evaluated in eq. |5]): 



'^P II cvap - ^/ 
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PfT 



A. 



GH 



(32) 



ai represents the contribution from the recoil of evap- 

II cvap 

orated particles to the variance of the longitudinal momen- 
tum; Aq-h has to be taken as in equations [TH] and [201 

Our final analytical formula is the quadratic sum of 
the different contributions, the Fermi motion (from eq. [21] 
for masses beyond ^lim or eq. [211 for masses below Aum), 
Coulomb repulsion between the multifragmentation prod- 
ucts (eg. 1261 applied below Aiim), and the recoil induced by 
evaporation (formula [321 for the whole mass range) . The 
longitudinal momentum dispersion of final fragments is ex- 
pressed by: 



P||Formi P||cvap 
''pilFcrmi ~'~ '''Pllcvap ~^ '^P||Coul 



beyond A\\ 



(33) 



below A\\ 



lim 



In fig. [21 we confront our predictions with three sets of 
data, measured in the reactions ^'^^Xe+Pb at 1 A GeV (56| . 
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Figure 3: Model in its full form as the square-root of eq. 1331 (plain 
line) and Morrissey's systematics (dashed line) compared to different 
data sets (23.[53.[5l. 

86Kr+Be at 0.5 A GeV [H and ^^Fe+p at 1 A GeV [El 
respectively. These data present longitudinal momentum 
dispersion of fragmentation residues over a broad mass 
range. Other published data concentrating only on very 
peripheral collisions (i.e. masses very close to the mass of 
the projectile) are easily reproduced by our model as well 
as by Morrissey's systematics. Since the present model 
does not consider the fission process, data containing fis- 
sion fragments were disregarded. 

The shape of the global prediction for cTpn of spectator 
fragments is in good accordance with the overall trend of 
the data available for this observable (see the confrontation 
with some data sets in fig. [3]). 

It is interesting to note that, although the abrasion pic- 
ture, which is a basis of the present model, is not expected 
to be valid for targets as light as protons, our model is 
able to reproduce such data (see fig..[3l) with a remarkably 
good agreement. 

The degree of sophistication seems sufBcient to repro- 
duce the data. Our analytical formula is in general accor- 
dance with the data, even for light fragments and no sys- 



tematic deviation is observed. In this mass region, where 
small residues can originate from more violent collisions, 
Coulomb interaction is reduced by the expansion of the 
system at break-up and our prediction could have been 
too large; this is not the case. 

The good agreement of our model with rather compre- 
hensive sets of experimental data give some confidence that 
the ingredients of our model may correctly represent the 
physics of the fragmentation reaction. A final conclusion, 
however, would require to investigate the sensitivity of the 
results to the ingredients of the model and to the values 
of the parameters. Such study is beyond the scope of the 
present work. 

5. Conclusion 

A new description of the width of the momentum distri- 
butions of fragmentation residues at relativistic energies is 
proposed. It describes the experimental data over a very 
broad mass range and includes most features of the known 
physical processes at play, constituting a great progress 
compared to previous formulae. 

The present new description has several advantages 
compared to predictions traditionally used. Morrissey pro- 
vided an analytical description of the data but only in the 
limit of small mass losses. Furthermore, his formula is em- 
pirical and gives no insight into the physics of the reaction. 
Goldhaber introduced the idea that the Fermi motion pro- 
vides the major contribution to the shape of the momen- 
tum distribution of the surviving fragments. The problem 
is that the quantity referred to in his model is not the final 
mass of the fragment which is observed. 

Going deeper into the description of the processes oc- 
curring in the collisions and using results provided by ex- 
periments, we suggested a new formula. It considers all 
relevant processes the system is undergoing in an analyti- 
cal formulation. This description reproduces the data very 
well. 

Our formula gives a good estimate of the kinematics of 
projectile-like fragments in the relativistic-energy domain, 
which is the energy range concerning the applications men- 
tioned in the introduction. Moreover, the physics basis 
of our description provides a realistic prediction of the 
longitudinal momentum width of fragmentation products, 
in the whole energy regime where the abrasion is a valid 
model for the collision phase, i.e. for energies above the 
Fermi energy. 
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Appendix 

To ease the use of the present model for practical applications, we propose in the present appendix the summary of 
needed equations, to calculate the momentum width of the fragmentation residues. We also list the needed parameters 
with their recommended values. 

The present model is valid over a large energy range (from about 50 A MeV to at least 10 A GeV), where fragmen- 
tation (and multifragmentation) processes occur. It applies to all kinds of projectiles (from light nuclei to uranium), 
independently from the target nucleus (therefore suited for any target, from protons to uranium). 



Input parameters: 

Ap : projectile mass 

Zp : projectile atomic number 

Af : final fragment mass 

The predictions are calculated for projcctilc-likc fragments in the projectile frame, but if Ap and Zp are replaced by the 
target mass At and nuclear charge Zt, the predictions apply to target-like fragments in the laboratory frame. 



Calculated output: 

fTpii ; dispersion of the linear momentum along the longitudinal axis for the final fragment (product of fragmentation 
reactions, excluding fission and fission-like processes). 



Model parameters: 

ro : radius = 1.4 fm 

: squared elementary charge = 1.44 MeV-fm 
u : atomic mass unit = 931.5 MeV/c^ 

AE : average induced excitation energy per abraded nucleon — 27 MeV 
< -Efermion > : mean energy of one nucleon from Fermi motion = 20 MeV 
K : inverse level-density parameter K = A/a = 11 MeV 
s : average consumed energy per evaporated particle ~ 15 MeV 

^kii"™ ■ average contribution of the thermal motion to the kinetic energy of evaporated neutrons or protons = 8 MeV 

Tbu : break-up temperature = 5 MeV 

Vexp = Vhu/Vo '■ volume expansion factor = 1.5 

Summary of equations: 

For a given projectile {Ap,Zp), we calculate: 

Pf = 281 • (1 - MeV/c , 



Ef = '^ MeV , 
2u 



V 



itherm 
kin 



_l_ / T?^ 



therm 
kin 




/ < Ek 



> . 



The mass limit for multifragmentation (see text) can also be defined: 

AE 



For Af > Aiim (final- fragment masses close to the projectile), multifragmentation did not occur. We need the following 
equations: 

eAj + AE ■ .4, 
= AE + e ' 



IGH 



2 ^2 



p% Agb.{Ap - Agh) 



Ap-l 



2 P%V^ ( 1 



At ^gh/ ' 
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The total momentum dispersion for this mass range is given by: 



For Af < Aiim, multifragmentation occurs. The following ingredients are required: 



Ephu - 2/3 ' 
•^exp 



Pf 

PFhu ^ 



•'exp 



12 \ -Epbu / 



_ K-AE 

K-AE + T^^ " ' 

Then we can calculate: 

'^P II Fermi 

And also: 

J2 



AbnJ 5 Ap-l V^GHy 5 Abu 



2 PfT ( 1 



Needed to calculate the contribution from Coulomb forces, the distance between the centre of mass of the produced 
fragment and the one of the whole system (i.e. the multifragmentation source) is taken to be 

''max = -Rbu — -RgH, 

with i?bu = ro (Au/Kxp)^^^ and -Rqh = ro (^gh/Kxp)^^^ ■ 
We define 

p2 0/1 AonA^^zy 2 (, Agh^' 
-n^oui = ^^GH • u -^j^3 Tjjjax • - -J— 

entering the expression of the Coulomb contribution to the momentum dispersion: 



.2 _ ^ p2 

■ -If 



1 

'Coul- 



PIICoul 3 

The total momentum dispersion for this mass range is given by: 

2 _ 2 . 2 I 2 

P|| ~ PlIFermi ~'~ P||ovap ' PUCoul' 
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